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ABSTRACT 

The GGSS is a partiaUy-fiUed, aU-sky survey to identify K-giant stars with low level of RV- variability. 
We study histograms of the radial velocity (RV) variability obtained in the early phase of the Grid 
Giant Star Survey (GGSS, Bizyaev et al., 2006). This part of the survey has been conducted with 
a very limited nubmer of observations per star, and rough accuracy. We apply the Monte-Carlo 
simulations to infer a fraction of the RV-stable stars in the sample. Our optimistic estimate is that 
20% of all considered K-giants have RV- variability under 30 ms^^. Different assumptions of intrinsic 
RV-variability for our stars give 12 - 20 % of RV-stable K-giants in the studied sample. 

Subject headings: stars: fundamental parameters, oscillations, late-type — techniques: radial velocities 



1. INTRODUCTION 

The future space interferometric missions have to be 
supplied with a reference frame of objects uniformly 
spaced in the sky, bright enough, and indicating no ra- 
dial velocity (RV hereafter) oscillations at a high level of 
accuracy (tens m/s). A large sample of relatively faint 
K-giants (V = 9.5 - 13.5 mag) that may serve as can- 
didates to the reference sources for the Space Interfer- 
ometry Mission (SIM Plan et Quest) has been studied 
in 2001-2004 bv IBizvaev eF a l. (2006 ) in the fram es of a 
preparatory s urvey GGSS (|P attcrso n~et al.ll200ib . As it 
was found bv IBizvaev et al . (2006), only about 1/3 of 
all candidates reveal the RV- variation within 100 ms~^ 
level, when observed with one year cadence. However, 
most of the objects in that sample were observed only 
2-3 times. Therefore, the real, or intrinsic variability of 
K-giants can not be estimated directly because of the 
limited number of RV-observations per star, and also be- 
cause of non-zero accuracy of individual RV estimations. 
A more accurate study of the RV -variation o f brigh t 
K-giants which was carried out by iFrink et ahl (|2001h , 
suggests a higher fraction of RV-stable stars among K- 
giants. On the other hand, well representative number 
of observations per star (from 11 to 28) rev eals 34 of 179 
RV-st able bright K -giants in the stud y by 
(12006!) ■ Studies bv lFrink et all (|200lD and 
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(12006) suggest a correlation between the RV variability 
and color (B -V). Contrary to them, this correlation was 
not found bv lBizvaev et al.l (2006), that can be explained 
by worse accuracy of RV estimates in the latter case. 

In this paper we undertake a study of histograms of 
the RV-variability p ublished for the GGSS survey by 
IBizvaev et al.l ()2006t ) . Monte-Carlo simulations that take 
into account the limited number of observations per star 
and rough level of the individual RV-estimations are per- 
formed. The main idea of the simulations is that the 
intrinsic (i.e. real and unknown) distribution of the stel- 
lar RV-variability is the same fo r both Northern (N) and 
Southern (S) subsamples (see iBizyaev et al.l I2006L for 
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detail). Since these two subsamples were observed with 
different instruments, and since different techniques of 
the RV estimation were applied, the observational un- 
certainties of the RV estimations in our simulations may 
be different for the subsamples. The subsamples differ 
from each other also in the number of observations per 
star. 

Resulting model histograms of the RV-variability (or 
sigma, as it was designated bv IBizvaev et al.l l2006l ) de- 
pend on both intrinsic RV-variability and factors intro- 
duced by observations. The published distributions of 
the RV-variability may be approximated by the model 
histograms by adjustment of model parameters. 

2. MODELING OF THE INTRINSIC SIGMA 
DISTRIBUTIONS 

We use the following basic assumpti ons to model 
the observational histograms presented bv lBizvaev et all 
([2006). 

(1) Distribution of the intrinsic RV-variability (or in- 
trinsic sigma, hereafter) is assumed to be the same for 
both N- and S-subsamples in our modeling. The shape of 
the intrinsic sigma distribution to be chosen later in this 
section. We note that distribution of the RV-variability 
is unknown for any limited number of observations. 

(2) Obsevations introduce uncertainties into the re- 
sulted RVs. We add normally distributed uncertainties 
to each individual RV. 

(3) We consider the reahstic sample sizes: Ngtars — 
148 and 341 for N- and S-subsamples, rspectively. The 
realistic numbers of observations per star Noi,s = 2 and 
3 for the N- and S-subsamples, respectively, are assumed 
in the simulations. 

For a certain value of the intrinsic RV variability, we 
generate Nobs normally distributed random radial veloc- 
ities (mean value of the radial velocity is zero, the stan- 
dard deviation is the intrinsic sigma) . Then we add a nor- 
mally distributed random " observational" uncertainty to 
each RV. We assume the standard deviation of this uncer- 
tainty as a free parameter in the modeling. The obtained 
model values of RV are used to calculate the mean ra- 
dial velocity and its standard deviation (i.e. the model 
sigma) exactly in the manner done with the real obser- 
vations (see the procedure described by IBizvaev et al.l 
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I2OO6D . The described procedure is repeated Nstars times. 
It gives us Nstars values of the observed RV-variabihty. 

The described simulation is repeated 40 times with ran- 
dom values of RV. For the comparison purposes, certain 
models are evaluated 100 times to make sure that re- 
sults are independent of the number of simulations. The 
resulting values of sigma are calculated for Nstars sim- 
ulations. The final RVs are obtained by averaging over 
all 40 repeated calculations. The obtained model sigmas 
are sorted out into the same bins in the same way as 
in his tograms obtained from observations (jBizvaev et al.l 
120060 . As a result, we obtain the mean numbers of model 
sigma in each bin of N- and S- histograms. 

For the comparison of the model histograms and ob- 
served RV- variability distributions, we figure out one chi- 
square value which is common for both N- and S- 
subsamples. The contributions of the samples into the 
are weighted proportional to the size of the samples. 

We consider the following three assumptions on the 
functional form of the intrinsic RV variability (sigma) in 
the modeling. 

(1) All stars have the same single value of the sigma. 

(2) There are two samples of stars: "RV-stable" and 
"RV- unstable", and hence there are two values of the 
sigma. 

(3) The intrinsic RV-variability is a continuous value 
which uniformly spans a range between and a certain 
value. 

The accuracy of the RV estimations introduced by the 
observations is a free parameter in all three cases. In 
the first case, an additional free model parameter is the 
intrinsic sigma. More flexibility is implied in the second 
case, where the "stable" (typically less than 50 ms~^) 
and "unstable" (typically greater than 100 ms^^) sig- 
mas, and fractions of stable and unstable stars are ad- 
ditional free parameters in the modeling. In the third 
case, the upper limit of the continuous sigma (the lower 
one is set at 0) is the free parameter. In the latter case, 
we approximate the continuous sigma by its 100 discrete 
values. 

3. RESULTS OF THE MODELING 

The modeling shows that agreement between the 
model and observing histograms of distribution of the 
RV-variability can be achieved with quite wide ranges of 
free parameters, whereas the fiting of the southern his- 
togram is more difficult. Our model #1 does not agree 
well with the observations. The best-fit intrinsic RV- 
variability in the model #1 takes high values of the or- 
der of several hundreds m/s, and the best-fit accuracy 
of both N- and S-observations is rough (several hundred 
m/s). The minimum value is about 3.5 times higher 
than that obtained in the next two models. 

The best- fit intrinsic sigma of the "stable" stars in the 
model #2 ranges from 10 to 30 ms~^, whereas the "un- 
stable" stars have the sigma of 890 ms~^. The best-fit 
accuracy of the RV estimations in this case is about 100 
and 50 ms^^ for the southern and northern samples, re- 
spectively. The best-fit fraction of the RV-stable stars 
is 20% in our model #2. Figure [T] shows the histogram 
based on observations (solid) and the best-fit model dis- 
tribution of the RV-variability (dashed). A higher frac- 
tion of the RV-stable stars would make the central peaks 
of the model histograms too high in both top and bot- 



tom panels. A lower fraction of the RV-unstable stars 
would create a disagreement in the tails of both N- and 
S-histograms. 

Model #3 gives the best value of all considered 
models. The best-fit observational RV-uncertainty is 50 
ms~^ for both N- and S-samples. The best-fit intrinsic 
sigma is distributed uniformly between and 850 ms~^. 
Assuming 100 ms~^ as the limit between the RV-stable 
and unstable stars, one can obtain 12% as a fraction 
of the RV-stable stars in our sample. Figure [5] shows 
the best-fit model and observing histograms of the RV- 
variability for the model #3. 

Although formally the minimum x^ value in the model 
#3 is less that that in the model #2, the difference is not 
significant. At the same time, the comparison of Figures 
[1] and [2] shows that the model #2 better matches the 
shape of the observing N- and S- histograms. 

The real distribution of the intrinsic RV-variability can 
have an arbitrary shape. More precise studies of the RV- 
variability suggest a one-peak distribution of the intrin- 
sic RV-variability and a higher fraction of the RV-stable 
stars than that found in our simulations. In order to 
make the fraction of the RV-stable stars higher in our 
modeling, one had to assume that the maximum of the 
intrinsic sigma distribution would have been shifted to- 
ward the lowest values of the RV-variability. On the one 
hand, our model #3 is an extreme case for such class of 
the intrinsic sigma distributions, and hence 12% is the 
lower limit for the fraction of stable stars. On the other 
hand, the model #2 is another extreme case, since the 
uniform distribution of the intrinsic sigma is a marginal 
case of arbitrary monotonously descending distributions. 
As a result, the real fractions of the RV-stable and un- 
stable stars would not be that much different from 20% 
and 80%, respectively. 

Any deviations from our assumptions about the model 
parameters may change the assessed fraction of the RV- 
stable stars. Thus, introduction of systematic errors into 
measurements of RVs makes the model histograms wider, 
allowing the fraction of bona fide RV-stable stars to be 
higher in this case. The observed RVs for the southern 
sample show worse agreement with our modeling than 
the northern RVs under any assumption of the balance 
between the fractions of RV-stable and unstable stars. In 
the case of the southern sample, we might encounter with 
systematic errors in the observing RVs. Note that the 
southern sample spectra were obtained w ith rather low 
S/N, of the order of 10 and even less (see iBizvaev et aP 
I2OO6I for details) . However, most of the northern stars 
were observed just two times. It introduces more uncer- 
tainties into the modeling for the northerm subsample. 

4. CONCLUSIONS 

Our Monte-Carlo simulations of observed distributions 
of the RV-variability reveals 12 - 20 % of the RV-stable 
stars in the 2001-2004 sample of the Grid Giant-Star 
Survey under basic assumptions on the intrinsic RV- 
variability distribution. Our optimistic estimation is 
that 20 % of all considered K-giants have the intrin- 
sic RV-variability under 30 ms~^. It is in agreement 
with a recent study of the RV-va r iabil ity performed for 
bright K-gian ts by Hckkcr et al. (2006), but less than 
that found by lFrink et al.. (2001 ). As it was discussed by 
IBizvaev et al.l ( 20061 ). constraining the physical parame- 
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Fig. 1. — Distributions of observed RV- variability (sigma) for the 
Northers (top) and Southern (bottom) samples are shown by the 
solid lines. The dashed lines connect the points on the model dis- 
tribution of the RV- variability. The latter distribution corresponds 
to the best-fit solution obtained within the model #2 (see text). 



ters of the considered candidates can help to increase the 
fraction of RV-stable stars, and partially explains the dif- 
ference in the stable stars fractions reported in different 
papers. 

The universal receipt for a radial velocity study in- 
tended to select the RV-stable stars is to conduct as 
many observations per star as possible, and to achieve 
the accuracy of individual RV- estimations that is bet- 
ter than the required level of the RV-stability. In this 
case, the RV-variability histograms obtained from the 
observations will reveal the shape of the intrinsic sigma 
distribution. It reflects the spirit of recent searches for 
the RV-stable stars for prosp ective space interferometric 
missions (Hek ker et ai]|2006( ), and agrees with the ideas 
and design of the ongoing SIM Planet Quest grid can- 
didates observations where the iodine cell technique is 
being applied. 
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Fig. 2. — The same as in Figure [Tl for the model #3. 
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